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TORNADOES AND OTHER ATMOSPHERIC VORTICES 

by RoIjcTt Cl. IX'issler 
LewiH Research Center 

Abstract 

'Hie ^rawth of rarKiom vortices in an atmosphei'e with buoyant instability and 
vertical wind shear is studied analytically. A study is also made of the veloci- 
ties in a single gravity-driven vortex; a frictionless adiabatic model which is 
supported by laboratoi'y experiments is first considered. 'I'he effects of axial 
drag, heat transfer, and precipitatlon-induccxl downdrafts are then calculated. 

Heat transfer and axial drag tend to have stabilizing effects; they reduce the down- 
drafts or updrafts due to buoyancy. It is found that downdrafts of tomadic magni- 
tude might occur in negatlvely-buoyant columns. 'iTie radial-inflow velocity re- 
quired to maintain a given maximum tangential velocity in a tornado is detei^ 
mined by using a tui-bulent vortex model. Conditions under which radial-inflow 
velocities become sufficiently lai-ge to produce tangential velocities of tomadic 
magnitude are determined. The radial velocities in the outer regions, as well as 
the tangential velocities in the inner regions may be large enough to cause damage. 
I'he surface boundary laye»*, which is a region where large radial inflows can 
occur, 18 studied, and the thickness of the radial-inflow friction layer is esti- 
mated. Finally, a tornado model which involves a rotating parent cloud, as well 
as buoyancy and precipitation effects, is discussed. 

1 Introduction 

(Joservations of tornadoes have been made for many years, the first known 
photograph having been taken in 1884 (Science, 1972; Fig. 1). The destructive 
aspects of tornadoes are well known. Still, there appears to be little agreement 
about their dynamics. The lack of agreement is. of courae. at least partly due 
to the devastating nature of tornadoes and the consi>quent difficulty of making 


meaningful measurements. Less devastating atmospherle vortlees sueh as dust 
devils and waterspouts might be profitably studied, but there Is no assurance that 
their mechanisms are the same as those of tornadoes. Thus theories and labora- 
tory experiments should play an Important role In tornado research. 

Capabilities for observing the stronger vx)rtlccs have, of course, been stead- 
ily Increasing (Preprints of Papers, Ninth Conference on Severe Local Storms. 
1975). This Is forturuite, since theories and lalx)ratory experiments by them- 
selves are Incufficicnt. Thus we may perform an amlysis of vortex behavior, 
devise an experimeni to test the analysis, and get agreement between the two. 
However, there is no assurance that the results will apply to real atmospheric 
vortices, unless we know that the assuniptions made in the theory and e\p* rim(.*nt 
apply to vortices In the atmosphere. In order to make progress, theory, lalK'ra- 
loi’y experiments, and field observations are all indispensiblc. 

Several aspects of atmospheric vortices are considered herein, including 
their production, structure, and maintenance. Although the emplvisis will be on 
toiTiadoes, much of the discussion will also be applicable to smaller atmospheric 
vortices, such as dust devils and waterspouts. 

Throughout the paper we will use the Boussinesq approximation in the cqua- 
tloas of motion of the fluid; that is, the fluid will be considered as incompress- 
ible, except where density differences affect the buoyancy, and the density differ- 
ences will be t;iken to depend only on the temperature. However, the effects of 
moisture and vertical pressure differences on the buoyancy can usually be accoun- 
ted for by using suitable (X'tential temperatures in place of ordiruiry temperalures 
in the equations (e.g. , Kuo, 19()6). 'Hie potential temperature is dcfini<l as the 
temperature attained by the fluid when it is compressed adiabaticaliy t< a stamlai-d 
pressure. The fluid may contain condensing or evaporating liquid. If the potential 
temperature is independent of altitude, there will be no buoyancy force acting on a 
displaced fluid element. Thus for the results to be applicable to atmospheric 
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vortices, the temperatures in the analysis should be considered to be potential i 

temperatures. 

In the toiTiado models considered, axial symmetr>’ will be assumed. AlthouKh 
there may be significant nonsymmetiic aspects of ioi-m»does, most of the important 
features should show up in an axisymmetric model. 

2. Growth t)f Random Vortices In An Unstable Atmosphere With V'ertiial Shear 

Extensive studies of tornado and funnel-cloud occurrences (Wills, 1969; and 
Modahl and Gray, 1971) have shown a high degree of correlation of those occur- 
rences with atmospheric instability and vortical wind shear. Wills has stated 
that it does not appear coincidental that the major tornado-producing region of 
the world (the eastern two-thirds of the United States) is the only area where 
strong vertical wind shears and strong potential instability are frc^iucntly simul- 
taneously present. ITius one approach to the study of tornado occurrences might 
be to analyze the growth of random vortices under these conditions. 

Such an analysis has already been carried out in a study of the growth due to 
buoyancy of a homogeneous tuiHaulence (or rortex distribution) with shear 
(Deissler, 1967. 1971). The self- interaction between the turbulent eddies (ran- 
dom vortices) was assumed to be small compart'd with their interaction with the 
mean temperature and velocit>’ gradients. (This assumption docs not necessanly 
imply that the random vortices are weak.) That case is of interest here, since we * 

are more concerned witli the interaction of the rortices with the mean gradients 
than we arc in their interaction with one another. 'Hie initial distribution of 
vortices is assumed isotropic, but their distribution becomes anisotropic under 
the influence of the mean gradients. Although this approach considers the oecur- 
rence of atmospheric vortices as random events, it should show conditions under 
which such occurrences arc highly probable. From a forecasting standpoint this 
may be all that cap be hoped for; the prediction with certainty of the exact time 
and place of a tornado occuricnce may not be feasible. 
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Hesults for the growth of the vorticlty anti size of typical vortices are shown 
In Fig. 2 for a Klchaixlson number HI based on vertical temperature and velocity 
gradients of -10 That Is a typical Hlchaixlson number for tornado occurrence 
(Wills, 1069), where (Hi a /3g(dT/dz)/(dU/dz)^) . (The vorticlty of a typical vor- 
tex is taken here as the itx)t mean square vorticlty at the wave number where 
Uie vorticlty spectrum Is a maximum. The size of a typical vortex is taken as 
the reciprocal of the wave number at which the vorticlty spectrum is a maximum.) 
'llie qiuintlty ^ Is the thermal expansion coefficient, g is the gravitational body 
force per unit mass, dT/dz is the vertical gradient of potential temperature 
(negative for unstable conditions), dU/dz is the vertical gradient of horizontal 
velocity (vertical wind shear), and Is a characteristic wave number (recip- 
rocal of size) for tlie initial vortices. 'Ilie value of xj shown in the figure Is 
somewhat arbltraiy but should be of the light oi-der of magnitude for an unstable 
atmosphere if the vortices under consideration aro themselves tuHiulent. In that 
case the kinematic viscosity u will be replaced by a much larger turbulent vis- 
cosity (Delssler and Peiimutter, 1960). 

The results In Fig. 2 show that both the vorticlty and size (particularly the 
vorticlty) of a typical random vortex Increase considerably with dimensionless 
time (dl'/dy » 0.003 sec”S. 'Hiis result is illustrated schematically by the ror- 
tlces sketched in Fig. 3 for an early and a later time. Thus, in agreement with 
tlu* observational studies of Wills and of Modahl and Cray, Uie theoretical results 
show that the simultaneous presence of strong potential instability and vertical 
wind shear arc favorable for the development of atmospheric vortices. The re- 
sults also show the importance of the time element in the strengthening of atmo- 
spheric vortices. In addition to the Hichaixlson number, a dimensionless time 
(dl7dz)t is also an important p:<rameter. (Alternatively we could, by multiply- 
ing (dl /dz)t by (Hich.aixl.son number)^^\ consider the dimensionless time to be 
( dK(dT/dz)^^“tJ. ) 'ITius if the conditions of strong iwtcntial instability and 


veitii'a! wind shear persist for a looij enough time, the oeeurrenec of strong at- 
mospheric vortiees will become highly probable. 

One point which might be emphasized is that in this model we have not had 
to postulate a large supply of cyclonic or anticyclonic vorticity. 'lire net voiliclty 
can in fact be zero, some of the vtjrticcs being cyclonic and some of them anti- 
cyclonic. On a large scale, this generation of rortices may be similar to the 
split of a thunderstorm into anticyclonic and cyclonic storms (Fu)ita and 
(irandoso, lUGW). Of course the large vortices (tornadoes or tornado cyclones) 
are more likely to be cyclonic because of the effect of earth rotation. Thus 
lines of tornadoes (cyclonically rotating clouds) tend to move away from squall 
lines (nonrotating clouds)(Fuglta, 1975). Uut the smaller dust devils are nearly 
as often anticyclonic as cyclonic. (According to this point of view the movement 
of cars on the left instead of the right side of the road, to generate less cyclonic 
and more anticyclonic vorticity (Isaacs et al, 1975), would not appreciably in- 
fluence the total number of tornadoes, but it may decrease the proportion of 
cyclonic to anticyclonic ones. In order to do that it may be necessary for the 
road to be hotter than the surroundings (Manton, 197G).) 

Can we conclude, then, that vortices of tornadic intensity can be generated 
and maintaintH.1 by the mechanism considerc'd in this section? The c ' t be 
generated in that way, but if they are like tui-bulent eddies, they will decay in 
times on the order of a few multiples of r/v, where r is the vortex radius and 
V is its velocity. Thus the Dallas tornado (Hoecker, 19U0) had a radius at the 
maximum tengential velocity of alx)ut 70 metera anti a maximum tangential veloc- 
ity of about 70 meters per second so that its lifetime, if it is like a tuii)ulent 
eddy, would be on the order of a few seconds, a time much loo short to be 
realistic for tornadoes. On the other hand, a vortex having the same value of 
radius times tangential velocity as the Dallas toraado (a measure of the vortex 
strength), but a diameter of several miles, could have a lifetime on the order of 
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u few hours. Thus the large to. 'yclone or rotating parent cloud (FuglUi, 
19tiU) from which the tornado contracted, might have been similar to a turbulent 
eddy formed in the way considere<i in this section. Hut in ortier for a vortex of 
tornadic strengUi and size to be mainUiincd. large rsidial inflows will be rc<|uircd 
(see section 5). 

Another point which should perhaps be discussetl is that buoyancy-irt' iced 
insUibility tends to produce hoHzontal vorticity (Deissler, 1967), whereas the 
vorticity in tornadoes tends to be more nearly vertical. However since only a 
small percentage of the vortices develop into tomuioes, it will be sufficient if 
at least some of them are n*;t horizontal, l erhaps the reason veilical shear, 
as well as buoyancy, Is apparently required for tornado development (Wills, 

1969) is that vertical shear tends to protluce more nonhorizontal vorticity than 
does buoyancy (Deissler, 1967). 

Once a vortex has been generated in an apparently random fashion, as dis- 
cussed in this section, we can consider its structure. We will now look at a 
vortex in more deUtil. 

:i. A Frictionless, Adiabatic Model 

For the sake of simplicity we will consider first a frictionless vortex with- 
out heat transfer. The analysis given here is a generalization of that given by 
Dc'issler and Boldman (1974). 'Fhe steady-state etjuatlons for the conservation 
of mass and momentum for an axially symmetric flow can be written, with th“ 
Houssinesq approximation, as (Lamiau and Lifshitz, 1959) 
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where r uml z are, respeetlvely, Uie radial and vertiral eoordinales, u, v, 
and w are the radial, tangential, and vertieal veloeltles, p is the density, p is 
the pressure. T is the potential temperature in the vortex, p^^ is the is|uilibrium 
pressutx', is the etjuillhrium temperature at z ® 0, j{ is the aeee'.enition of 

uravity, and /j » T ^ is the thermal expansion coefficient for the atmosphere. 

Ihe last term in Kt|. (2) is the buoyancy term, where the vertieal temjn riture 
(’.radient in the surrouiKlin^ atmosphere dT^/dz, and T, are assume<l eonsUint. 

Let w = w(z> and v “ v(r> within the vortex or reuion of interest. If the first two 
terms in K(|. (-ii are small eompaied with the third term, as they piobaljlv are 
for most of a tornado dloeeker, IDtin, Kq. (1) shows that p - p^, ip - p^,Mn. (In 
the outer region, althounh the first teim may not be small eompareil with the thirci, 
p - p^, as calculated fiom Ki|, (1) either with or without the first term, is small.) 
Kquations (2) and (.'b become 
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Solution of K(i. (5), with the initial condition w = w^^ when z = U then nives, 
for the vertieal velocity. 
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w - i ; w“ - - T)z - dU ^ z‘ 

dZ 
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From Kqs. (1) and (7) we get. by using the condition that u is finite at r = o. 
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u « T . 




- 2<in(T^.y - T)/. - ,JK(a i^/i»z)?.- 


11/2 


Kqiuilions (7> and (8) show that for larKc z, 


(8) 


W * ±(-/iK dT /dZ)^^“z; u» (-/fK dT /dZ)^^^ r 
e ^ c 

M 

whereas for small z, 

w « Wy-, u » */>K<Tpy - 0 r/(2Wy) (10) 

Thus the radial veloeity u ean be lai-Re If Wy is small and T^^ - T is non- 
zero at z ~ 0. As will be seen in section 5, a large radial inflow is favorable for 
the maintenance of large tangential velocities. 

From Kq. (G), we get 


v/vj » (r/rj)"^ (11) 

where Vj is the tangential velocity at r » I’j. This equation applies reasonably 
well except for very small r, where friction beromes imporUint (I)eissler and 
l•el■lmutter. lOGO). The important question ol when tangential friction cannot be 
negleetcfl will l)c considered in detail in section .'i, 

.Although the vertic-al temperature gradient d'r^/dz is important in the at- 
mosphere. i*8 ef'ect is hard to investigate in the labonitory. 'I'hus, in orcier to 
see if the model considered here can give reasonable correspondence with lalK>i*a- 
lorv results, we devised experiments where the buoyancy force was proportional 
to T^,^^ - T (or to - p) . For that case, along a streamline (Deissler and 
iiohiman, 15/74). 

u ^ dr _ 1 r 

w dz 4 7. 
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or 


- 1 J X 

r/Pj - (z/zj) (12) 

where z^ is the value of z at Home inner reference radius r^. A plot of Kq. (12) 
is shown in Fig. 4, w’herc Wq is 0 and z increases neKativcly downwartl. 'Hiis 
plot is not intended to imply that the velocities in a tornado are necessarily down- 
ward; the stream lines could just as well be plotted converfilnR upward. How- 
ever, in our experiments the vortex was heavier than the surroundinR fluid, and 
its motion wan downward. 

In one experiment (Ueisslerand Uoldman, 1974) the boiloff from liquid 
nitix}Ren, to which a small amount of initial vorticity was added, provided a source 
of cool, heavy gas which could accelerate downwanl to produce a concentration of 
vorticity Conden.sation streamers made the flow visible. 

Top and side views of the resulting gravity-driven vortex are shown in Fig. 5. 
'Ilie predicted stre .mline shown in Fig. 5(a) was calculated from (Deissler and 
Holdman, 1974) 




V’l 


(r/Py)"** - 1 
(r/ry)-* - 1 


(13) 


where the central angle \p was .set i>qiuil to zero at an outer radius Py and 
is the value of (p at an inner radius Pj. Both the theory and the experiment 
show a characteristic hook shape for the streamlines. In the outer portion of 
the vortex the flow is nearly radial, whereas in the central portion it beromes 
nearly tangential. 

It might be mentioned that the vorticity concentration shown in Fig. 5 was not 
dependent on having a flat surface immediately above the rortex. However a 
hollow cylinder several vortex diameters high and c-losed at the top was placed 
above the vortex to reduce the vertical velocity Wy to a low value at the plane 
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when' thi* fluid cnlcTfd the vortex. Aeeordlnn to Kq. (10), then, the radial veloi- 
ity at that plane waa lar)te, and vortieity ivneentration eould take plaee. In an 
aetual tornado a temperature inveraion (aUible reKion alxive the vortex may have 
somewhat the same effeet (Wa ixi.l 050) . Similar results wt . e obtainiHl when a 
plate was plaeed immediately above the vortex. The plate simulates the ground 
li the upward and downward direetions are reversed, .so that the flow is thought 
of us uein )4 upward. 

The side view of the vortex shown in I-'i);. 5(b) appears to be similar to that 
predieted in Flu. ■*. Jif least for the rejcion reasomibly near the top. However 
to show the vertieal developnu*nt of the vortex somewhat more elearly, and to 
investigate the effeet of vertical ‘5f„'ai in the enviionmcnt on the development of 
the vortex . we (Holdman and Deissler) devise*! anothi>r experiment, in which a 
robttin)’ column of dyed dark corn syrup (speiifie jjravity. 1 . .1) accelemted tkiwn- 
waixl thix)uj>h horizontally flow ing water, in this ease the vertical veloi-itv at the 
point where the heavy fluid entered the water was restric*<-d to a comparatively 
low value by the friction in a loUilin); tube through which the fluid flowixl befo»-e 
it enlercxi the water. 

Huns were made with and without lolalion and with and without water flow . 

The results for die four cases are shown in Fl^. G. (In order to show the columns 
more clearly, the corn syrup lyinn on die liottoin has been removed.) 'I'he con- 
traction of the column as it accelerates ilownward to produce vortieity concentra- 
tion (in the rotatin); ca.ses) is clearly .seen in this ex)K*riment, and is similar to 
flat (iredlcted in Hr. 4. The results with and without roUttion ajiiwar to be about 
the same (except that the columiis in the roflitin)t cases whi)) around more). This 
re.sult miRht be exjiected. if the analvtical model in this .section is applic.able, 

since Kq. (5) or (7) is indeiwndent of rotation, and since the continuity condition 

*} 

for this ca.se can be written as D*^ c consflint. Thus the diameter I) of the 
column at any vertical position should be the same with rotation as without. 
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C'ompu riHon of thf rt'tiultH with and without horl/.unUil wator How (vtM'tlcal wind 
Hhuar) Meama to ahow that the* columns with water flow arc somewhat smaller than 
those without. One mlKht be inclined to expect this intuitively, since he* night expect 
the vertical shear to cause additional stretching of the column. If, however, one 
nuMiifics Kq. (.*>) for the case where the column makes an angle 0 w ith the vci>> 
tical. Kq. (5) becomes, 

/ dT 

tk V dZ ' 


where w an z are in the axial direction, and Z la in the vertical direction. 
But since dZ » dz cos d, Kq. (14) becomes 


wilJi--)lg T - T i:z 
dZ dZ 


dT 


( 1 51 


which has the same fonii as the wuation for d - u if w and Z are. re.spcctlvcly. 

•> 

in the axial and vertical directions. Since the continuity condition is still l)“w’ -* 
consUint. the column diameter with and without horizontal water flow should be the 
same at a paiticular vertical position. Therefore an\ tl life rent- es between the cases 
with anfl without horizontal water flow in Fig. G. if Uiey arc rial, must be due to 
the effects of axial friction. Those effects were neglected in K(|. (15). 

Thus, at least for these e,\pei1ments. ami with the possible exi i*ption of U>e 
effects of friction, the model considered in this section seems to be realistic. 

For a vortex- to-cn\ii*onment density ratio closer to one than the value of 1.3 in 
the experiment, as it probably Is for a tornado, the axial drag may be relatively 
more imporUint. The effects of drag, as well as several other effects, will be 
considereil in the next section. 


( 
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4. Kffecis Of Drag. Heat Tranafer, And Precipitation- Induced IX^wn Drafts 

On 'Ilie Axial Motion 

We first estimate the axial drag on the vortex column. The axial turijulent 
shear stress on the column (molecular shear stress is negligible away from walls) 
is given by 

T,, ■ -p w'u' « -p (w - U sin «)" 

U I il ty 

(w - L’ sin 0)“ 

* -p_H(w - I' sin 0)^ sgn(w - U sin 0) (Hi) 

where p^ is the air density, w' and u' are turbulent velocity components in the 
axial and radial directions, the overbar designates an averaged quantity. 0 is the 
angle of inclination of the column with the vertical, w is the velocity of the col- 
umn relativ.' .i. .. stationary point, and U sin 0 is the component of the hori- 
zontal wiiiu velocity along the axis of the rolumn. 'I*hus the axial velocity of the 
column relative to the surrounding air is w - I’ sin 0. According to experiments 
of Liepman and Laufer for a free turbulent mixing layer (Townsend, 19.56), 

H iJillHlJ 0.01 (17) 

(w - I' sin 0)“ 

The symbol sgu in Kq. (16) means sign of, and is included so tluit will 
change sign if w - U sin 0 changes sign. 

If rain drops arc present in the column, they also exert a force if they move 
relative to the air in the column. The force of the drops acting on unit volume of 
air is 


F 


C - d* 

Lij n 


•? Pn" 2 

•” (w^i - w) 8gn(w^j - w) 


(IH) 


where Cjj is the drag coefficient of the drops in the air, d is the drop diameter. 
Wj is the axial velocity of the drops ri'lative to a stationary point, and n is the 
number of drops per unit volume. 
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If wo w hto an axial force balance on a volume element of fluid of diameter I) 
(ineluutnK buoyancy, ptvelpiUition, and axial dra^ forces) , where I) is the column 
diameter at axial position z, we cet. in place of Kq. (14), 



where N Is the number ol rlmps passinK thi-ounh un(t cross-sectionjil area of the 

column in unit time and is related to n in Kq. (18) by N » nwj. Z is measui-ed 

•> 

in the vertical direction, z in the axial direction. If drops arc conserved, I)“N 

is constant for steady state. It has been assumed that the vertical wind Knidient 

is constant, and that I’ * 0 for Z » 'Ilius I' has been ri*placed by ? 3L'/yZ 

•> 

in the last term in Kq. (19). Hy continuity , the quiintity D“w in Kq. (19) is a 
constant, since for realistic cases the cross sectional area of the drops Is small 
compared with that of the column, 

'Hie first (i‘:’m on the rinlu side of Kq. (19) nives the effect of buoyancy on the 
axial acceleration, the second term nives the effect of rain diops or hall, and the 
last tei-m five's the effect of turbulent drag. 

In order to get an equation for Wj, we write a force balance on a dion. This 
gives 


2 

dw <1 f- 1) f>,, •) 

= -2g cos 0 ~ (w . - w)“ 8gn(w . - w) 

dz 2 d '^d “ ” 


( 20 ) 


where is the density of the drop. 

If heat transfer takes place between the column and the surroundings, the 
temperature T in the buoyancy term in Kq. (19) will vaiy with z. 'Khe turbulent 
heat iiansfer q^^ to the column per unit area can be related to the shear stress 
by a form of Hcynolds analogy (IXdssler, Ui-W): 


H 


c (T - T)1 t,,| 

pe 0 '*0' 

I w - (dl’/dZ) 2 sin 0 \ 


( 21 ) 


where Cp^. is the effective specific heal in l ie column, and the quantity in the 
denominator is a((ain the axial velocity of the column at z relative to the surround- 
inj{ air. (We do not have to evaluate Cp^,. since it drops out of the final equation.) 
Writing an energy balance on an elemcn*. of fluid of diameter D, and using Eqs. 

(16) and (17) for the shear stress in Eq. (21), give 


- T) 
dz 



w - — Z sin 0 
dZ 


( 22 ) 


Before calculations can be made, we still need a relation for 0, the local 
angle of the column with the vertical. In oixler to obtain such a relation, we 
assume that U - u' « constxint, where u' is the hoiizonUil velocity of the column 
at z. so that du'/dZ = dU/dZ. Using that relation and the relation u' ■ w sin 0, 
we get 


dz dZ \ 2 / dz 


( 2 . 1 ) 


2 2 

from which 0 can be obtained as 0 * (w 6)/w . The horizontal and vertical co- 
oi'dlnatc.'^ (X, Z) of the column at z arc then given by 


— = sin 0, — * cos 0 
dz dz 


( 24 ) 


2 * 

Tl»e diameter D of the column at z can be obtained from D = (D“w)/w, where 
D“w is a c*onstant of the motion. 

We can now solve the set of ordimiry differential Eqs. (19). (20), and (22) to 
(24) numerically, if we know the initial conditions anti the constants in the equations. 
For the initial conditions we set w, w^. - T, w^O, X, and Z equal to 0 at 

z = 0. For the constants let (i = 1/.130 g = 9.81 m/sec^, dT /dZ = -0.00.1 K/m, 


V... 
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C|) “ 0. 4 (for a solid sphere!, d* 0.001 m, H«0.01, dl'/dZ 0. 00.3 see” ^ , and 
* 0.0008G. The poteniial temperatures and wind veloeltles eorrespond to 
average values glvxm by Wills (1900) for tormido environments. 

Caleulated results for D“w * -10** and -10^ are plottinl In Figs. 7 ami 8. 

■> U 

For eomparlson, we note that the Dallas tornado hud a D~w % 4x10 (lloeeker, 

1960). Aeeoixling to the results of Wills (1969), the effeets of buoyancy and vem 

tleal wind shear In tornado- producing environments extend over the vertical dls- 
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lanees shown. Values of U“N of 10 and 10 In Figs. 7(a) anti (b), respectively, 
with drop size d » 0.001 m, eorrespond to violent rain with .some hail (Ileri^ 
et al. , 1945). However, results calculated foi no rain agreed within a percent 
with those In Fig. 7. Thus the effect of rain in producing dowmlrafts in a tornado 
appears to be small eomparetl with the effect of buoyai;oy for those uilues. How- 
ever preeipltatlon-indueed downdrafts may be important for eonveeting lai*ge-seale 
vxirtleity downward at a slow rate in cases where buoyancy effects are small 
(Fugita, 1975). Calculations for no buoyancy but with the precipitation parameters 
the same as In Fig. 7 gave air-downdraft velocities of 0.5 m/sec witli a drop veloc- 
ity of alwut 7 m/sec at Z « -6,000m. (In oixler to calculate nonzero velocities from 
Kqs. (19) and (20) for this ease it was necessaiy to assume small initial uilues for 
w and w^|.) Wlien the drop size was increased an order of magnitude to 0.01 m, 
a size which might occur for hail, downdraft velocities approached those for nega- 
tive buoyancy (I'ig. 8). Giant hail is in fact observed in connection with tornado 
occurrences (Browning, 1965). For (his case it was even possible to obtain 
downdrafts with small positive buoyancy (dT^/dz = ^d.000:i). 

The contraction of the column as it is accelera -* downward by buoyancy forces 
is clearly shown in Fig. 7, and is similar to that in the experiment in Fig. 6. 'Hie 

relative contraction for a given change in altitude is greater for the larger value 
o 

of D“ w because the effects of axial drag and heat transfer are less for the larger 
column. As is the case for the drag, the effect of heat transfer is to decrease the 
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vertical acceleration. ITils latter decrease occurs because the buoyancy tenn In 
Kq. (19) Is re<luced by huU transfer to the column. ITie vertical wlmJ shear de- 
flects the column horizontally as It falls. It appears that the main effect of wind 
shear on the axial acceleration of the column Is to decrease the axial drag term 
(last teiTO in Kq. (19)). In no cases calculati‘d was the effect of wind shear large 
enough to reverse the sign of the axial drag term. 

Values of axial velocity are plotted against vertical position Z In Fig. 8. 

'fhe rasults extend over a wider range of values of D^w than those in Klg. 7. 

The smallest value of D^w may be appropriate for dust devils as well as for 

9 

small tornadoes. The axial velocities decrease as D“w decreases because the 

restricting effects of drag and heat transfer are greater for the smaller eolumns, 

as discussed in the last paragraph. 'ITie curve for D“w » « is the same as that 

for no drag and hejrt transfer. Comparison of the solid with the dashed curve 
6 

for D“v/ » 4x10 shows that the vertical wind shear increases the axial velocity 
at a given vertical position, 'fhls increase Is due to the drag reduction produced 
by vertical shear (discussed In the prcctxling paragraph). 

For the lai'gest columns, the vert. cal velocities near the ground (If the 
ground is at Z =-8,000 m) approach tornadic values on the oixler of 70 m/sec. 
Kven higher values would be obtained if the dc>wndraft originjited at higher alti- 
tudes, If the local potential temperature gradient were negatively greater Uiun 
-0.003 (as it may well be), or if giant hail wero present in large quantities. Thus 
aside from the effect of tangential velocities (to be investigated in the next section), 
it appears that buoyancy- induced downdrafts may be lai*ge enough to pioduce sig- 
nificant damage. This effect has been discussed by Kossman (1960). lie presents 
a photograph by K. L. V'an Tassel (U.S. Weather Bureau), shown in Fig. 9, which 
gives evidence of severe downdrafts in a tornado which passed thmugh Bayaixl, 
.Nebraska in 1951. Pieces of 2- by 8-lnch lx»ards were driven into the ground a 
depth of 18 inches. The ends driven into the ground were .square (not splintered). 
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Fufclta (1975) has also noted evidence of severe downdrafts in tornado environ- 
ments. These observuvions appear to be consistent with the ealeuiated results In 
Hg. 8 for the larger columns. 'Hius at least some of the damage produced by 
tornadoes appears to be due to downdrafts. 

If we neglect the preclpltation-drag term and let - T*0 at z»0 in 
Kq. (19), there is nothing in that equatiou which says that the motion should be 
either up or down; the buoyancy force in that case is zero at z = 0. Actually, 
from the observations and analyses of various authors, it appeal's that the veloci- 
ties in some tornadoes (or parts of toi'nadocs) may be up, whereas in others it 
may be down (Hoosman, I9(i0; Fugita, 1973, 1975; Davles-Jones et al. , 1975; 

6 

and Golden, 1975). Since Fig. 7 was for downflow, an upflow case (D“w = 4xio ) 
is plotted in Fig. 10. Most of the discussion for Fig. 7 applies also to this case. 

As pointed out by Kskridge and Diis (1970), and by Dclsslerand lloldman 
(1974). one difficulty with models which consider only updrafts is that the toi-natk) 
would be expected to form before the rotating cloud forms, whereas the opposite 
is observ’ed. As will be seen, there are probably regions of Ixith updrafts and 
downdrafts in tornadoes. 

5. Maintenance Of The Tangential Velocity 
Except for a frictionless vortex model, the tangential velocity has not yet 
been considered. The effect of turiiulcnt friction is important, in that it determines 
the maximum tangential velocity which can occur in a vortex; the maximum tangen- 
tial velocity is often used as an indication of the violence of a tornado. Obviously 
the frictionless relation v « cannot apply all the way to the vortex center. 

In order to investigate the tangential velocity, wc use Burgers' (1948) model 
for a steady-state visc-ous vx>rtcx, as modifletl for a tuHjulcnt viscosity by Ueissler 
and Berlmutter (19G0). The modifie<l etjuation for tlic tangential velocity is 
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where and Vy are respectively the radial and tangential velocities at an 
outer vortex radius r^, and k is the KdrmAn constant. For most of the cases 
considered here, the flow is radially Inwaixl (Uq Is negative). Hy continuity, 
when a radial flow is present there must also be an axial flow. In deriving 
Kq. (25), we let w var>' with z only (w « z), and u/Uy * *'/•*(>,• These relations 
(Including Fq. (25)) satisfy the equations of motion and continuity, 'rhe outer vor- 
tex radius I'y can be considered as the radius where w in an actual vortex be- 
gins to fall off rapidly with increasing r. ''.*he molecular kinematic viscosity 
was replaced by a turbulent viscosity t , and t was obtained from a modified 
form of von K.irmfln's similarity hypothesis as (Deissler and Perlmutter, 1960) 

e = «^Vyry/2 (26) 

where v was taken proportional to r \ as it is for the outer portion of the vor- 
tex. The constant value for t given by Kq. (26) was assumed to apply to the 
whole vortex, since the turtjulence (at least for most of the cases of interest) is 
convected from the outer regions to the center of the vortex by the radial inflow. 
Strictly speaking, Kq. (25) applies only to the case where u and v are iixle- 
pendent of z. However, outside of the regions of steep gradients in a boundar\’ 
layer, Kq. (25) should be sufficiency accurate for the purposes for which it is 
used here. In Deissler and Perlmutter it was found that for confined vortex 
flows, the KdrmAn constant k»! 0..'1. That value of k was also found to be 
reasonable for the experiments of Wan and Chang (1972) for a tornado-like 
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vortex In the ivKion above the Hurfaee boundary layer. 'ITius the value of x « O.H 
Is retained here. 

A plot of Eq. (25) Is ){iven in Fig. 11. In the outer region of the vortex, va r”^. 
whereas for small r, v * r. As -Uy Increases, the region for v oc extends 
to smaller r, so that the maximum tangential velocity increases (for a given Vy). 
Thus for large tangential velocities, large radial inflows are required, fora com- 
parison a case where the radial flow is outward (Uq/ V y“ 1) is also shown. Tine tan- 
gential velocities are small except at large radii. This kind of velocity distribution 
may occur in rotating clouds, where the flow may be outw'ard (Ilsu and Fattahi, 15175). 

We want to relate the maximum tangential velocity to the radial-inflow 
velocity Uj^ at the same point. To that end set dv/dr from fk). (25) equal to 
7^ro. linen let Uy » u^ Tq/ r^^, since (according to our model) u varies linearly 
with r. and let r * r^^, the radius at which v is a maximum. Tills gives 
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E(iuat ion (25) becomes, for r = i- and u, , = u r„/i* , 
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1 - exp'd/x Mv^r^/vyry)(u^/v^)(ry/r^)- 


For large (greater than about 2.5. which will occur if v « r 

near I'y) these two equations give, with x = 0..1. 
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Thi8 rema I'kably aimple relation allows us to estimate the radial-inflow veloeity 
required to maintain a nlven maximum tanKt'ntial veloeity. ’I'hus for the 
Dallas tornado is ^iven by Hoeeker (I960) as about 70 m/see at a radius r^^ of 
about 70 m. These values, or values slightly higher, probably afiply also to the 
Cleveland tornado (Lewis and Perkins, 1953). Thus from Kq. (27) u^^ * -10.8 
m/sec. 

We now ask whether a buoyancy- induced sink is large enough to produce this 
required va’ue of radial inflow velocity For no axial friction and no heat trans- 
fer to the column (these factors would tend to decrease u^), and for T^,y - T •» 0 
at z » 0, and Wq *= 0 (or for large z), Kq. (9) gives, for Will's tornado- 

environment ciata and for the I3allas tornado, u^ « 0.33 m/sec. This is smaller 

m 

than the above value of u^ required to produce the observed tingential velocity 
by a factor of ulx)ut 33. 'Ilius it appears that buoyanc-y effec'ts far from any re- 
striction of vertical-flow velocity (large z) are much too small to sustain the 
tangential velocities observed in typical tornadoes. Acc^ording to Kqs. (9) and (27) 
the potential temperature gradient would have to be about (33)“ times as great as 
the average value observed for tornado occurrences. However for veiy large 
vortices, with r^ on the o»xler of a mile or more (large for a torreulo), Kqs, 

(9) and (27) show that tangential velocities of to roadie magnitude might l)c main- 
tained by buoyancy forces in regions away from vertic-al-flow res trie-lions. 

(Jnc possibility we should investigate is that the large tangential \elocltles 
in tornadoes may be produced transiently, say by random turbulence effects. 

(In a random distribution some eddies will be very strong.) In that case we can 
estimate the oitier of magnitude of the dec-ay time t^j (with small radial flow) 
from the unsteady equation for the tangential velocity (Goldstein, 1938) as 

V V 2 V 

iv _!n V r 

ait, 2 2 2 
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where the kinematic viscosity was rc*plact*ci by a turbulent visc-osity t Klven by 
Kq. (26). For the I>allas tornado t^ ^ 20 sec, a time which, although longer than 
that obtained earlier by simply considering the vortex as a turbulent eddy, is still 
too short for a toi*nado. (Average tornado lifetime is about 13 minutes (Kskridge 
and Das, 1976).) 

Thus it appears that the Uingential velocity in tornadoes must be maintained 
by radial-inflow velocities on the order of those given by Kq. (27). Uut the val- 
ues given by Kq. (9) for regions away from vertical-flow restrictions, are much 
too small. Kquation (10) shows, however, that for regions where Wy is small, 
as above a surface (for upflow). or beneath a stable inversion (for downflow), 
u can be large if - T is nonzero at z = 0. Thus we let Wy = 0 and 
T^ - T = 0, -2 and -1 K, and r = r * 7»> m in Kq. (8). 

cu ni 

The resulting curves for u^^ are plotted against z in Fig. 12. The large 
effect of T^y - T near z = 0 is evident. At least for a few meters aloove the 
surface (or beneath an inversion), the values of u^ are of the right oi*der of 
magnitude to maintain the tangential velocities observed in tornadoes, according 
to Kq. (27). Thus a vortex potential temperature which is a few degrees different 
from the surroundings in a region where the vertical velocities arc restricted to 
low values appears to be necessary- for maintaining the vortex by buoyancy forces. 
Conditions particularly favorable for this type of vortex maintenance should occur 
for fire whirls (and possibly for dust devils) where high local surface tempera- 
tures are prevalent (high T^^y - Ty). Of course, in those cases, the smallness 
of r in Kq. (8) would tend to woi'k against niaintenance by buoyancy. 

This section has emphasized the tangential velocity (maintained by the radial 
inflow) as the destructive element in tornadoes. According to Kq. (27) the radial 


vrloclty iH only about 15 poivent of iht* UinKintlul velocity ut the nulluH where the 
UinKontial velocity Is a maximum. However, tlie i-atio - (u^ 

varies from 0.5 to 1..T for the ralues of u^/vy In Fig. 11. Thus the nidial velocity 
In the outer part of the wi’tex can be of the same oixler of maKnitude as the maxi- 
mum tanitential velocity. 'Hits appears to be the case for the experiment in 
FIr. 5a, where the streamllnc.'i are nearly radial in the outer part of the voHex. 

Accoixlinfi to F(|. (8). u becomes infinitely larRe at a surface where z » 0. 

Wy “ 0, and T^,y - T / 0. In rcsility, of course, u decreases in a thin boumUiry 
layer alxjve the ground and goes to zero at the ground. The important subject of 
the surface Iwundary layer will be considered in the next section. As will be 
seen, the surface boundary layer can have a large effect on the radial inflow, even 
when buoyancy effects near the surface are small. 

0. The Surface Boundary Layer 

'iTie surface boumlar>’ layer below a vortex has been studicxl analytically by 
Kuo (1971). Lewd len (1962). and others, and expcrimenUilly by, for Instance, 
Diang (1969), Wan and Chang (1972). and Savino and Keshock (1965), According 
to the physical picture usually given, the centrifugal force produci‘d by the ro- 
tating flow Is in approximate c<|uilibrium with the radial pressure giadicnt in the 
region above the Ixiundary layer. Within the Ixiundary layer the centrifugal force 
is reduced because v is smaller there. If the pressure is approximately con- 
stant across the boiimlary layer, so that the radial pressure gradient does not 
change appreciably, there will be a net inward force in the Ix'undary layer which 
lends to produce lai*ge radial inflow's there fhese radial inflows (given by Kq. 

(27)) may be large enough to mriintair the tangential velocities obsi*r- i-d In tor- 
nadoes. The intensification of the vortex b\ the pri'sence of the boundaiy layer 
has been shown by Hsu and TalUihi (1975), who ran i-xpcriments both with aral with- 
out a ground plate (no buoyancy). Waterspouts, which travel over water, may not 
be as violent as tornadoes, which travel over land, because the tangential velocity 
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probably doe:> not go to zero at the surface of the water, as it does at a land 
surface. 

After the fluid leaves the boundary layer vertically, it will continue upward 
if the local buoyancy force is large enough (see preceding section). Otherwise 
at least part of it will tend to recirculate, as in the experiments of Wan and 
Chang (1972) (Fig. 13). It is of inteiest that in those experiments the airflow 
in the outer part of the vortex was downward, even though the air entered the 
apparatus radially from a location below the ground plate and exhausted through 
the ceiling. 'Unis recirculation will no doubt take place in a toiTiado unless the 
potential temperature in the vortex at the ground is sevei'al degrees higher than 
that in the surroundings, as in Fig. 12. 'Hiere is some evidence that air 
descends from a rotating parent cloud at large ladii (as a collar cloud) and then 
travels radially inward in the surface boundary layer (e.g. , the Xenia tornado 
(Fugita, 1975)). Fsk ridge and Dus (1976) have indicatated that precipitation- 
induced downdrafts surrounding tornadoc's may be warm. In that case the air at 
the surface which enters the tornado updraft may have a potential temperature 
several degrees higher than that of the surroundings, so that according to Fig. 12 
and Kq. (8), radial inflows sufficiently laiyie to produce UingerUal velocities of 
tornadic magnitude (Kq. (27)) may be productni by buoyancy forces. Although 
this is an appealing suggestion, its occurrence will di*pend on the presence of posi- 
tive buoyancy, and on precipitation drops of sufficient size and density to bring 
down the wai-m air. A sample calculation using Kqs. (19) to (2-1) showed that with 
the conditions for the precipitation case in Fig. 8 (dot-dashed curve), but with 
positive buoyancy (dT/dZ = +0.0003), air which descended 8,000 m to a surface 
because of falling precipitation was about 1.3 K warmer than the surroundings 
at the surface. Its downdraft velocity near the surface was about 40 m/sec. If, 
on the other hand, the air descending at large radii is cool, it will have to be 
pushed up mechanically at smaller radii. 
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In c‘aHt'8 where the buoyancy foreea are Hmull, what is the source of enerK>- 
for the tornado? Apparently it is the rolatinK parent cloud almve the v’ortex. 

Thus the rolatinK cloud centiifugr's air outwani which then descends, swirlinK 
at larKc radii. At the surface at least tuirt of the air becomes trappinl in the 
bound:iry layer, whei'e it spirals inward, and then upwani, to form the vortex 
(FuKlla, 11)75; ’*».u anti Fattahi, lt)75; ilsu, 1972; and S/illinsky and Wippcrman, 
1969). (This is somewhat similar to the stirred teacup effect.) Fxperimcntal 
flow patterns obtained by Hsu and Fattahi (1975) for a lotatinK honeycomb above 
a surface are shown in FIk. 14. .Note that the vortex pjotiuced at the surface is 
off-center; it follows a circular path aix>un<l the centerline of the apparatus. 

There is also recent evidence (e.K. . the Xenia tornado (FukIUi, 1975)) that 
considerable damage is tione by small, intense short livisi vortices, called suc- 
tion vortices, which drop to the ground from the lati;er vortex. 'ITie surface 
layer probably also has an Important effect on the dyruimics of these smaller 
vortices. 

Next we estimate the thickness 6 of the friction layer where u goes from 
zero at the ground to its maximum value. Alwve that point u again decreases. 
Calculated values for 6. using tangential velocities and radii for tornadoes, 
should at least show whether such conditions are likely to occur for a boundary 
layer. 

The efjuation for the boumlary layer for the radial flow is 


•> 



dr i)z r p dr pr dz 


(29) 


where p is the pressure, t is the shea r stress, and Vj^ is a value of v in the 
frfetion layer. At some distance outside the friction layer for u the centrifugal 
force temi will t>e in approximate cquilibrtum with the piessure-gradient term, or 

iiJ2 

r p d r 


( 30 ) 


25 


An assumption K*-*ne rally made in analyses of vortex bountiar)' layers Is that the 
pressure Is lonstant from the wall to the dlstanee where Kt|. (3U) applies. U'e will 
retain that assumption, althouKh the data of Wan and C'hanK (1H72) indieate that It 
is only approximately eonstant. Also, the thlekness of the frietion layer where u 
Koes from zero to its maximum x^tlue is small eompared with the thiekness of the 
boundary layer for v (Wan and C'hanK. 11^72), so that « u in Kq. (2ih. Kqua- 
tlon (29) beeomes, with K’q. (30), 

u . w ♦ X (31) 

Or Oz r pr Oz 


where v is mi'asured outside the boundary layer, where Kq. (30^ apnlies. 'Ibe 

first term on the rlKht side of Kq. (31) ace«'len>tes the fluid radially inward in the 

frietion layer. In the remaln(k‘r of the analysis we i-onslder quantitieK at r^^. 

the radius where v outside the Ixrundary layer is a maximum, At that 

radius ou/or u_/ r„ , Moreover w varies linearly with z for small z. 
m m 

Also, sinee the boundary' layer is tuHjulent, the veloeity u parallel to the sur- 
face can be approximated by -u « . Hy in*iodueinK these approximations 

into Kq. (31) and inteKratinK from the wall to z where -u is a maximum and 

T “ 0, we Kft 
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where t^., the shear stress at the wall, can be approximated by the Hlasslus 
relation 


— = 0.0228 u^ 
p m 



(33) 


where e is the kinematic visioslty. By performinK the same operations on the 
continuity Kq. (1). we get 



'Hu- ratio *** nu-ato ml at t - wiuTfa.< that in K«|. (l!7» 

i^ nu-aHuml near Uu- ulKt* of Uic boumlary layi*r for v. Slmx* Ihr valui* of 
calculated from Fq. (.'15> 1 h somewhat InscnMitlve to it shouhl be* suffi- 

ciently accurate for our puipose to use the value from Fq. (27>. Settinu 

-u / V * 0. IA5 iilvcs 
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Note that the surface i)oundai'> layer is Uu* onlv part of a toi-njulo where the 
kinematic viscosity r is nut ncKliKil>lc compared with the turbuU*nt viseositv. 

For Case 2 In Wan an<l (’hand's (iy72> experiment (Case 2 has features simi- 
lar to those of a tui-nadoK their Tal)lc 1 t>ives r^^^ “ U,07U m and v^^ * l!».M 
m/sec. With i- ■ I . .!■ 10“'^ m“/ sec, Fcj. '"!i> <;ives ‘*1 1 . a value that 

appears consistent with results in Wan arsl Chant's Fin. !-• Similarly for the 

Dallas tormdu. wheie r » 70 m and v ~ 70 m/sec (Hoeiker lOGO) we net 

mm 

= 0, IG m. nils value, althounh small compared with the otlier lennths In the 
tornado is, accordinn to K(|. (MG), larne enounh to maintain the Uinnential vtducity 
by radial inflow Althounh -u ({uiekly reac-hes a maximum as z liu-rejises to 
its flecrease with further increases of z Is, aceoniinn to tlie experiment of 
Wan and ('hann, much more nradual. ■Iluis, tlie radial velocities In the lx)und:»ry 
layer should be hi^h etwunh to mainLxin tai.nential velocities of tornadie manni- 
tude for .some dlsOincc al)Ove the ni''>und. Fejuation t.'t.'i) is for a smooth surface. 
Fora rounh terrain 5^^ would be larner than the value nlvin by Kq. (.'U>) 
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7. Stitnmary And ConclusionH 

AlthouKh we have not been able to obtain a aingie model applicable to all torna- 
doea it may be that there la. in fact, no unique mechanlam for all tomadoea. We 
have pointed out probable waya the strong vortices may oi1glna»e and persist, as 
well as how they may be structured. 'ITiesc suggestions were supported, as much 
ns possible, by calculations and observations. 

Ihindom vortices were shown to grow with time in an atmosphere with buoyant 
instability and v’ertical wind shear. 'Fhus if those comiltions persist long enough, 
the occurrence of strong atmospheric vortices becomes highly probable. 'Hie 
analysis showed that strong radial Inflows must occur in the vortices, if their high 
tangential velocities are to be maintained for a reasonable length of time. 

A relation between U»e maximum tangential velocity in a toriiado and the 
radial-inflow velocity at the same radius (Kq. (27)) was obtained from a tuiHaulent- 
vortex model. The radial velocity was found to be rough’; 15 percent of the tan- 
gential velocity at the radius where the tangential velocity is a maximum. It was 
shown that tangential velocities of tornadic magnitude could be maintained by 
buoyancy-induced ladial inflows. However, it would be necessary for the potential 
temperature in the vortex to differ from that of the surroundings by several de- 
grees at a point where the vertical velocity is small or zero (at the ground surface 
or below a stable Inversion). Buoyancy effects in regions away from surfaces or 
stable layers wero too small to maintain large tangential velocities, except for 
very large vortices with radii for the maximum tangential velocity on the order of 
a mile or more (large for a icruudo). 

If buoyancy effects are small, the boundary layer next to the ground (where 
centilfugal effect are rcxluced) must be relied on to allow the large radial inflo.. s 
rcHjuired to maintain the vortex. The thickness of the friction layer where the 
radial velocity goes from zero to a maximum as height increases, was found to be 
extremely thin 'less than a meter) for a smooth surface. For a rough terrain the 
thickness would be greater. 
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A physl'^iil picture which seems to emei’Kc. and which may be appiicable to 
many tornadoes is the following A larj;e \orte\ or I'otatinn cloud (turl)ulent edily) 
is produced, possibly by the effects of buoyancy and vertical shear, as discussed 
in section 2. The rotation of tlie parent cloud tends to centrifuge air in the out- 
waixi direction. 'Hie rotating air then descends at large i-adii oecause of the 
vertical circulation pattern set up by the rotating cloud. 'Hits descent may be 
aided by negative buoyancy effeet« and precipitation, if the potential tempc •ature 
at the cloud is lower than that at Uie surface. If, on the other hand, the potential 
temperature at the cloud is higher than that at the ground, buoyancy will tend to 
rotard the descent at large ra'’ I bat the air may still descend because of ;.iv- 
cipitation and the natural vertical circulation pattern set up by the rotating cloud. 

As the rotating air reaches the ground, at least part of it will be trapped in 
the surface boundary layer and will flow inward at high radial velocities (high 
enough to produce tangential velocities of tornadic magnitude), 'iliis radial in- 
flow is caused by the reduced centrifugal force in the lx)undai7 layer and, if the 
air is warmer than the surroundings (if it came from a .'arm cloud), will be 
aided by the buoyancy sink produced by rising air at small radii. If, on the other 
hand, the air is cooler than the surroundings, (he radial inflow will be hindercxl. 
'Hie air will flow upward at small radii because of the mitural circulation |:>attern 
set up by the rotating cloud and, if the uir is warm, because ot buoyancy torccs. 
Of course if the air is cool, its upwaixl How will be hindered. As Uie air reaches 
the cloud level, the cycle may be repeated. 'Ilu* tornado would b«' e.xpected to last 
as long as the rotation of the cloud continu<*s to be sufficiently large and/or buoy- 
ancy and precipitation forces continue to lu- effective. 

According to this model a b)nuido mignt be pr'jduced regardless of which way 
the buoyancy force's act. The essr ntial ingredient seems to be the roUilion of the 
cloud. Precipitation at large radii slu)uld also have a positive effect regardless 
of the dirc'ction of buoyancy loives. The' small inti'iise suction vortices ebser' ed 
in tornadoes may also have essentially the mechanism outline :l here. 


AlthounH muc-h of tht- da' :.uki> done by tornadoes is ptx>i>ably euused by their 
hinb tangential veloeities, ealeitlalions ami observatioas show that downdraft veloei- 
ties jn nejtatively buoyant eolumns ean reaeh toniadie ma(>niludes. The analytieal 
model, which was supported by laboratory experiments, took into ueeount the 
effects of buoyancy, axial drUK, heat transfer, and precipitation flow. Both axial 
di-aK ami heat transfer tended to reduce the downdraft and u(Miraft velocities 
si^tnificantly . (They had stabili/Tn(> effects.) Wlien buoyancy was present, pre- 
eipitation liad a comparatively net^li^ible effect wlien values of dro|) flux and size 
for violent rain were uswi. However when Ijuovancy effects were absent, the pre- 
cipitation causeti tile column to still descend at a low velocity. Wlicn dix)p sizes 
an order of ma^cnitude hiRher were used (1 cm), as mi);ht occur for hail, the down- 
draft velocities pixniuced by precipihition flow appmaih those produced by i)uoy- 
ancy, 

Tiie results also showt^l that i-adlal-inflow velocities in the outer part of a 
vortex may be compai'al)le to Uie maximum tangential velocities. Allhoujtti the 
radial inflows and downdrafts in a tornado may cause considerable destruction. 

U.e destruction associated with the tan);ential velocities may be still j>reater, !'(•- 
cause of the low central pressuros which accompany the tangential velocities. 
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Figure I. - First tornado cloud knorin tohase tieen pnotografiiiMt iMiner 
County. South Dakolo. August. 1 884. 
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Figure I • Groarth of vorttces in an unstaMe 
atmospnere with vertical shear. Heavier 
lines indicate stranger vortices. 


Figure < - Predicted variations, anth vertical distance. o( radius 
and tangential velocity along a streamline lor a tnctionlcss 
adiabatic model. 
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(at Column not rolitinq. No horitontal water flow. 


IM Column not rotating. Water flowing horizontally 1?$. 4 cm-sK*^. 



(cl Column rotating (tube atxive column rotating (d) Column rotating (lube above columti rotating at 2.3 rpil. Water flowing hori' 
at;.3rpii. No horizontal water flow. zontally(25.4cm-»ec'l|. 

Figure 6. - Photographs showing (townward acceleration of a column of heavy liquid (specific gravity • 1. 31 In water. 
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f Igurt 7. • Plots showing contrKtIon ol column 
KC*l«r«ting downward with nogjtiv* Ixioiitncy. 
vortical wind snaar, and prKipitatipn. ifT/dZ • 
-aOOJ°k-m'\ dUidZ • d003sfc‘*. initial 
T, • T • a Drop sue d • D 001 m. Diamotar 0 
IS to same scale as other distances. 
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VERTICAL DISTANCE Z. Km 

F Iggre B * Variation al aaial velocity with vertical position, 
0^, potential temperature gradient, vertical wind shear, 
and prKipitation. Other quantities same as In figure 7. 
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fiQure q. - Evidence ot severe downdraft (photograph taken l»y f . 1. Van 
Tassel, ‘I. S, Aeather Burejui, 


VERTICAL »CIGHT z. 



CCXUMN RADIUS D/2. Km 

Flqurt Id • Plot showtnq cwtrKtlon ol column Kcoltftlinq 
upwofd «itn buoyoncy. D‘ w • 4x10^, dU/02 * d No pro - 
cIpiMion. Oinor quontitioi torn# ot In ligurt 1 . 
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Mqurt IL • Virwtion ol Unqinliol vtlocily with radius and 
radMl volKily. Subscripts 0 indicatt valuts al an outer 
raitrtnct radius 
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Fqurt 11 * Plot Showing rtcirculation currants in eipari- 
mtnt ol Man and Chang (197?) lusa 2i. r„ is radius 
where tangential yelocity is a maiimum. 
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RADIAL VELOCITY u„. m-SEC* 

Figure 12 - Variation ol radiai-inilow velocity u„ <at 
r^. where v is a maiimumi with vertical height and 
temperature difference af i * d as calculated from 
Eg. Ibi. wg • 0, dTj'di • -d 003, r^ • 70 m. 


aNKRUNt 

4 



f iqure U - (>p(nn>«nUI IIcm iMiterns 
obUinMl by Hsu *n<) f tttihi (197^1 lor 
i routing hontytomb ib»t i surlKi 


